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Abstract 

The polar condensation/sublimation of CO2, that involve about one fourth of the 
atmosphere mass, is the major Martian climatic cycle. Early observations in visible 
and thermal infrared have shown that the sublimation of the Seasonal South Polar 
Cap (SSPC) is not symmetric around the geographic South Pole. 

Here we use observations by OMEGA/Mars Express in the near-infrared to detect 
unambiguously the presence of CO2 at the surface, and to estimate albedo. Second, 
we estimate the sublimation of CO2 released in the atmosphere and show that there 
is a two-step process. From Ls=180° to 220°, the sublimation is nearly symmetric 
with a slight advantage for the cryptic region. After Ls=220° the anti-cryptic region 
sublimation is stronger. Those two phases are not balanced such that there is 22%±9 
more mass the anti-cryptic region, arguing for more snow precipitation. We compare 
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those results with the MOLA height measurements. Finally we discuss implications 
for the Martian atmosphere about general circulation and gas tracers, e.g. Ar. 
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non-condensable gas 



1 Introduction 



The "cryptic region" is a dark region covered by ice in the South Polar Region of 



Mars (Kieffer et al. (2000); Titus et al. (2008)). It appears to be a region where 



most of the "spiders" features are located (Piqueux et al. (2003); Kieffer (2007)). In 



order to simplify, we define in this article the cryptic region from longitude 50°E to 
230°E (through longitude 90°E), and the anti-cryptic region from longitude 130°W 
to 50°E (the complementary sector, passing through longitude 0°E) (see fig. [I]). 

During the phase of CO2 frost accumulation, the Seasonal South Polar Cap (SSPC) 



is mainly formed by direct condensation, but some snow events can occur (Forget 



et al. , 1998). GCM studies show that a topographic forcing by the Hellas basin 



creates an asymmetry in the mode of deposition. Precipitation events are more fre- 



quent for the anti-cryptic sector than for the cryptic sector (Colaprete et al. 2005 



Giuranna et al. |2008 ). This suggests that the texture should be more granular 
(smaller grain sizes) for the former compared to the latter. This texture difference 
produces a relatively higher albedo in the anti-cryptic region and also possibly a 
higher accumulated mass. 

The direct measurement of the sublimating CO2 mass on Mars has been done 



only recently by three different techniques: gravity (Smith et al. (2001); Karatekin 



et al. (2006)), neutron flux (Litvak et al. (2006)) and gamma ray flux (Kelly et al. 
(2006)). But the time and space resolutions of these methods prohibit any conclu- 
sion at the regional scale. At the present time, only direct measurements of the 
seasonal variations in altitude by MOLA are able to estimate the local CO2 mass 



sublimation (Smith et al. (2001); Aharonson et al. (2004); Jian and Ip (2009)). 



But this impressive technique is limited by low signal to noise ratio because sea- 
sonal topographic variations are typically less that one meter! To enhance the 



climatic signal, Aharonson et al. (2004) use a Fourier transform and apply a filter 



on the MOLA time variation keeping only the annual period. Results show that 
the amplitude of this annual Fourier cycle is asymmetric, with a maximum in the 
anti-cryptic region (see fig. 4 (a) in Aharonson et al. (2004)). This results lead to 
two different conclusions: (i) asymmetry in amplitude is actually related to asym- 
metry in accumulation and/or (ii) the annual period filtering may be not relevant 
to describe the real annual mass variation, so that the amplitude asymmetry may 
be an artefact related to the time of presence of CO2. The second hypothesis is 
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Figure 1. Definition of the cryptic and anti-cryptic region within this article. 



supported by the fact that the sublimation is much faster in the cryptic region 



(Kieffer et al. (2000); Langevin et al. (2007)) leading to lower amplitude of the 



annual Fourier cycle. If this hypothesis is valid, the bi-annual cycle should be 



stronger for the cryptic region. More recently, Jian and Ip (2009) focus on the 
south polar region and show that cryptic and non-cryptic regions did not present 
any significant thickness difference. 



Indirect measurement of the CO2 mass has been performed for decades through 
modeling of: (i) the recession of the surface covered by ice observed by visible/infrared 



detectors ( Leighton and Murray 


(1966) 


Paige and Ingersoll ( 1985 ); Lindner 


(1993); 


Forget 


(1998); Kieffer et al. 


(2000); Schmidt et al. ( 


2009c)), (ii) fitting 


the annual 


surface pressure cycle ( 


Pollack et al. 


(1990 


); 


Wood and Paige (1992 


); 


Hourdin 



et al. (1995)). The first modeling methodology has the advantage to allow simple 
estimation of CO2 ice sublimation mass at regional scale. It assumes that the CO2 
ice is in vapour-pressure equilibrium with the local atmosphere, and that the CO2 
mass balance at the surface is controlled by radiative balance. We propose to use 
this methodology using the latest OMEGA observation to estimate the local CO2 
mass balance in order to derive the difference in the cryptic/anti-cryptic region. 
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Figure 2. Definition of inner crocus line, outer crocus line and snowdrop distance. Stat- 
istics scale length L s t a t define the super pixel resolution, binning several OMEGA pixels. 
Black line represents the actual surface proportion covered by CO2 ice for each super 
pixel. Doted lines represents the arctan model law (Kieffer et al. ( 2000| )). Crocus line is 
defined in the inflection point of the arctan law. Inner and Outer crocus lines are defined 
by the surface proportion of 99% and 1%. 



For each considered location on the map, the so-called crocus date is the date 



(Kieffer et al. (2000)) when the surface deposit is completely gone. The crocus 
date is then a geophysical field of importance for atmospheric studies. Due to the 



fact that all instruments have a limited spatial resolution, Schmidt et al. (2009c) 



has introduced both inner and outer crocus date. Each location is defined by a 
super pixel with a spatial extent denned by the L sta t scale length (with a the pixel 
scale length L p i X <L sta t)- The inner crocus date - respectively outer crocus date - 
is the date when the first - resp. the last - pixels fails to exhibit any ice signature 
within the considered super pixel (see fig. [2j). 



We will use the D-frost model ( |Schmidt et al.| (|2009c| [2008J)) to estimates the C0 2 
frost mass balance. Several inputs are required: inner and outer crocus dates, al- 
bedo, altitude, slope. This model is able to calculate the total mass sublimated, 
from the date it becomes unstable (determined by the model itself), to the crocus 
date (determined by the observation). Thus the initial mass, before net sublima- 
tion, is not required as an input but will be determined by the model. 



The aim of this work is to: 



estimate the main input fields, required for the D-frost model, from OMEGA 
observations: albedo, inner crocus/outer crocus date, 
calculate the sublimated mass of the SSPC using D-frost. 
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• derive the height variations, assuming a constant density, and compare to MOLA 
measurements. 

• estimate the spatial asymmetry of the CO2 mass release. 

• derive potential effects on circulation at global scale, and on atmospheric gas 
tracers. 

To avoid confusion about the sign of the "sublimation rate" of the SSPC, we adopt 
two natural conventions: surface convention with negative mass rates (sublimation 
of seasonal frosts); atmosphere convention, with positive mass rates (pump up of 
the atmosphere). 



2 Methods 



We propose to use the OMEGA imaging spectrometer (Bibring et al. (2004)) on 
board Mars Express for simultaneously detect CO2 ice on Mars and estimate the 
albedo under the aerosols layer. A compilation of all dataset available is out of the 
scope of this article but should be an interesting step for the future. 



We use the following notation: latitude A , longitude 6 and time t. The spatial 
fields (inner and outer crocus dates) are sampled in a spatial grid of 0.3° latitude 
(200 steps from A = —90° to -30°) and 10° longitude (36 steps from 9 = 0° to 
360°). The spatio-temporal fields (albedo and sublimated mass) are sampled in 
the same spatial grid and the temporal grid of 5°Ls (72 steps from t = 0° to 
360°). Each element of the grid is called "bin". The succession of bins in latitude at 
constant longitude is called a "longitude sector". Respectively a succession of bins 
in longitude at constant latitude is called a "latitude sector" 

Due to the lack of the spatial/temporal cover we need to interpolate the measured 
crocus dates and albedo. The first two sections present our strategy. Each of these 
interpolation methods will be evaluated simply using the standard deviation of the 
difference between the actual data and the interpolated field. 



2. 1 Crocus dates field 



For each bin with available OMEGA observation, we estimate the inner and outer 



crocus dates using previously described method in Schmidt et al. (2009c). This 



method uses WAVANGLET, an automatic detection algorithm (Schmidt et al. 
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(2007)), and a special algorithm to extract crocus lines. We use the OMEGA 
dataset in year 2005-2006 (Martian year 27) from Ls = 110° to Ls=320° and 
selected 544 observations from latitude 30° S to 90° S. 



We propose here the following methodology to get an interpolated map of both 
inner tj n (A, 9) and outer t out (X, 9) crocus date. 



2.1.1 Interpolation method 

The crocus dates, binned into our grid (see sec. [2|, are re-estimated by a weighted 
sum of actually measured neighbouring values in order to interpolate the field in 
case of missing values. The weights depend on the distance and are chosen to be 
in a Gaussian shape. This operation is performed through a convolution product 
on the complete latitude/longitude grid with the Gaussian kernel: 

tZ terp (X,9) = ff t in (\9).- 1 .exp f^^£ + ^f^) d\'.d9'(l) 

JJ ^■O'lat-O'long \ ^lat ^i ong ) 



With : 

• computation on the complete grid : latitude x longitude bins: 200 x 36 

• variance of the Gaussian in latitude x longitude bins: o~i at = 41.66 latitude bin 
(12.5° latitude), ai ong = 5 longitude bin (50° latitude). The ratio of latitude and 
longitude variance is set to real average distance in our grid. 



2.1.2 Results and evaluation 

The interpolated crocus line field t™ terp , average for the cryptic and anti-cryptic 
regions are drawn in fig. [3] (a). The recession is clearly symmetric before Ls=140°, 
latitude 55°S. This symmetry can be biased by the interpolation, especially due to 



Hellas and Argyre basins where persistent CO2 ice have been detected by Giuranna 



et al. ( 2007a[ ). After Ls = 190°, latitude 62°S, the cryptic region crocus lines are 



20°Ls earlier than the anti-cryptic region. Two possible interpretations will be 
discussed in section|3j lower accumulated CO2 ice mass and/or lower albedo leading 
to higher absorption and sublimation. 

To evaluate the quality of this interpolation, we analyze the difference between 
actual measurement and interpolated field. The standard deviation is quite large 
(14°Ls for outer crocus date and 17°Ls for inner crocus date) due to the sparseness 
and discontinuity of the dataset. Also it could be due to the sensitivity of daily 
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variation in frost cover, ignored within our scheme. The mean snowdrop time 
is 15°Ls in agreement with the value measured by consecutive observations by 
Schmidt et aT] Q2009c[). 



2.2 Albedo field 



Radiative balance requires the directional-hemispheric bolometric albedo of the 
surface A, or simply "albedo", that is an integrated quantity of the BRDF (Bid- 
irectional Reflectance Density Function). It is dependent of the incidence direction 
(the solar direction) but integrated on all emergence directions of the top hemi- 
sphere and spectrally integrated - weighted by the solar spectrum. Definitions can 



be found in Hapke , 1993 



This quantity can be estimated from space using the reflectance spectra in the 
visible/near infrared recorded by OMEGA, by removing the dust aerosols contri- 
bution that significantly modifies the apparent surface albedo, notably for bright 



icy surfaces that appear darker. Vincendon et al. (2008) propose a method to eval- 



uate the contribution of the aerosol scattering and attenuation at all wavelengths 
from a narrow spectral region where the entire signal comes from the aerosols 
(saturated bands of CO2 ice). 



We use 265 images of the Martian year 27 from Ls=120° to Ls=300° in the latitude 
range 52°S to 90°S. We average all estimations points where the method provided 



by Vincendon et al. (2008) can be applied to the corresponding bin. Then we 



perform an empirical spectral integration of the solar spectrum in the visible and 



near-IR as described in Schmidt et al. (2009c). Finally, we estimate the complete 



field of seasonal deposit albedo using the following empirical interpolation method. 
The main advance of this interpolation method compared to our previous study 



(Schmidt et al. (2009c)) is to take into account both space and time evolution. 



2.2.1 Interpolation method 

The albedo, binned into our grid (see sec. [2]), are re-estimated by a weighted 
sum of actually measured neighbouring values in order to interpolate the field 
in case of missing values. The weights, depending on the distance in space and 
time, are chosen to be in a Gaussian shape. Albedo is interpolated in a similar 
strategy to crocus date (section 2.1.1) with a convolution on the complete lat- 
itude/longitude/time grid with a Gaussian kernel. In order to keep in a realist 
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situation, we apply a threshold on the resulting albedo and then smooth the res- 
ults. 



'1) Convolution with a Gaussian kernel on the complete latitude/longitude/time 
grid: 

A interpl (\,6) = {{I A(\, 9). ' 



J J J 27T .(Jlat-O long -(J time 

fix -A') 2 (e-e 1 ) 2 (t-t') 2 \ IW , , 

exp y — '- + K — - + K - dX'.dO'.dt' (2) 

\ ^lat 2,&long ^°~time ) 

With : 

• computation on the complete grid : latitude x longitude bins: 200 x 36 x 72 

• variance of the Gaussian in latitude x longitude x time bins: Oi at = 41.66 
latitude bin (12.5° latitude), ai ong = 5 longitude bin (50° latitude), atime = 4 
time bin (20° Ls). The ratio of latitude and longitude variance is set to the 
real average distance ratio in our grid. 

(2) Application of a threshold on the albedo A mterpl value in the range [0.3, 0.9] 
to get A interp2 to remove noisy data. 

(3) Convolution to smooth the threshold (same than eq. |2| using A mterp2 to get 
A interp , with: 

• computation on the local grid only to save computation time: 8x2x2 

• variance of the Gaussian in latitude x longitude x time bins: ai at = 8.3 
latitude bin (2.5° latitude), <Ji ong = 1 longitude bin (10° latitude), o~time — 1 
time bin (5° Ls). The ratio of latitude and longitude variance is set to the 
real average distance ratio in our grid. 



2.2.2 Results and evaluation 



The interpolated albedo field A mterp , average for the cryptic and anti-cryptic re- 
gions are shown in fig. [3~| (b). Albedo increases in both region of the SSPC as 
previously noted (Paige and Ingersoll (1985); Kieffer et al. (2000); Schmidt et al. 



(2009c)). The anti-cryptic region has an albedo 0.03 higher than the cryptic in the 
beginning of the recession. The difference increase to reach 0.2 around Ls=220° 
and then decrease again. The main reasons for a lower albedo in the cryptic region, 
compared to the anti-cryptic region should be the following: (i) the fact that the 



ice is translucent and thus the low albedo is due to the underneath ground (Kieffer 



et al. (2000)) (ii) the presence of the fans associated with the spiders (Piqueux 
et al.| Q2003|); |Kieffer et al.| (|2006[); [Kieffer | Q2007[)) (iii) dust contamination close to 



the surface (Langevin et al. (2006)) (iv) larger grain size in the case granular ice 



(Warren et al. (1990)). Recent studies suggest that a slab ice is present in extensive 
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(a) Crocus lines 



(b) Albedo 



Figure 3. Averaged fields interpolated from OMEGA observations for both cryptic and 
anti-cryptic regions, (a) crocus lines (b) albedo under the aerosols. The interpolated fields 
are used as an input of the D-frost model. 



position in the SSPC with a fast space and time evolution (Schmidt et al. (2009a)). 



We will discuss these points in an upcoming paper (Schmidt et al. (2009b)) 



This interpolation method is centred because the convolution with a Gaussian 
kernel is normalized. The standard deviation is 0.028 and is one order of magnitude 
lower than the mean interpolated albedo field (0.405). The resulting empirical 
interpolated albedo field is thus very close to the original measured dataset. 



2.3 D-frost model: daily average sublimated mass model 



We use the D-frost model to estimate the daily average CO2 mass surface balance 
dM = dM(X,8,t) for each bin 0.3° latitude, 10° longitude and 5°Ls as described 
in introduction of section |2j This method estimate the following radiative balance 
in the visible and IR range: 

dM = (F t t rm - WZn - WT cat - WT R - W}l or - F™ nd )-^—.dt (3) 



F°herm represents the thermal energy flux from the surface; F^ n is the direct solar 
energy absorbed by the surface; W™ at is the solar energy scattered by the atmo- 
sphere and absorbed by the surface; Wf^ represents the thermal energy emitted 
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by the atmosphere and absorbed by the surface; WJ? oor is the thermal energy emit- 
ted by the neighbouring surface facets and absorbed by the surface; F™ nd is the 
heat flux from the regolith; Lco 2 = 590 x 10 3 J. kg -1 is the latent heat for CO2 
sublimation. 



The thermal energy flux is estimated with the frost temperature determined by the 
local pressure of CO2, based on the LMD GCM (Forget et al. (1999)) - calibrated 
to Viking Lander pressure curves -, scaled to MOLA elevation. The albedo in 
the thermal infrared domain is assumed to be constant at 0.1 and the emissivity 
is set to 0.99. Despite these parameters are actually changing in space and time 
( Eluszkiewicz and Moncet ( |2003 )), a quantitative sensitivity test has shown that 
it should play a minor role (see section 4.1.6. of Schmidt et al. (2009c)). Also the 



annual heat wave is neglected to speed up the computation. This energy source has 
also been demonstrated to be negligible during the sublimation phase (see section 



4.1.3. of Schmidt et al. (2009c)). More details about the computation scheme can 
be found in Schmidt et al. (2009c). 



2.3.1 Ar description 

We improve the D-frost model by including the effect of the non-condensable gas 
on the CO2 partial pressure. We use the measurements of Ar by ISprague et al. 



(2007) to correct for the actual CO2 pressure using the formula: 

Pco 2 = Pgcm — Pat (4) 

In other words, we use Pco 2 instead of Pgcm in the estimation of surface temper- 



ature (see eq. 13 of Schmidt et al. (2009c)). 



We use the interpolated albedo field presented in the last section. Altitude and 
slopes are estimated for each bin averaging a map in stereographic south projection 
of MOLA data at a resolution of 920.8 m. 



2.3.2 Time integration 

The total mass sublimated M is estimated by this integral: 

-90° 360° t end 

M= [[I dM(X,e,t).dS.dt (5) 



-30° 0° t 
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where to is the beginning of the sublimation (i.e. the first time when dM < 0), t en d 
is the crocus date (i.e. the time when all the C0 2 ice disappears). 



We propose two extreme cases for t en d- the inner crocus date t{ n and outer the 
crocus line t out . A most realistic case is the decreasing of surface proportion X = 
X(t) inside each bin from the inner to the outer crocus date. Then eq. [5] becomes: 



-90° 360° t end 

M= J J J dM(X,e,t).X(t).dS.dt (6) 

-30° 0° t 



We assume a linear decrease of X(t) with time, from X = 100% at ti n to X = 0% 
at t out . Thus the surface proportion can be written: 

X(t) = (7) 

tout <^in 



In the following, this case will be referred to as "linear surface proportion". 



2.4 MOLA height variation 



The MOLA observations were performed from February 28, 1999 and June 30, 
2001 in a time interval longer than one Martian year (MY 23 and MY24). In this 
study, we used MOLA data from orbit 10012 to 20327 (from Ls = 103° to 360° 
in MY23 and then to 190° in MY24). Between Ls=344° in MY23 and Ls=16° in 
MY24 there are no MOLA data acquired because solar conjunction. This time 
interval permits to investigate the height changes of the CO2 polar deposits. 

To determine the thickness of the CO2 frost layer in different regions as a function 
of time, it is most important to compute the crossover residuals resulting from 
subtracting the altitude h of a certain location at some value of Ls to the reference 
altitude ho at the same location without frost. In this work about the south polar 
region, the no frost reference altitude is chosen at Ls=300°. We have limited our 
consideration to a latitudinal band between 50°S and 85°S. Several small areas 
of about 0.1°-0.2° in longitude and 0.05° in latitude and with enough number of 
crossovers are then chosen. After computation of the average elevation value of 
each small area, we calculate the residuals by subtracting the average altitude h 
of a certain location at some value of Ls to the reference altitude h$ at the same 
location. 
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3 Results 



The total sublimated mass from latitude 50°S to 90°S during the year is about 
5.9 x 10 15 kg (fig. [4]). The two extreme cases (4.3 x 10 15 kg for t end = t in and 
6.9 x 10 15 for t en d = t out ) are considered as the error bar due to the main uncertainty 
on the crocus dates such that our estimation is (5.6 ± 1.3) x 10 15 kg. 



Our estimation is compatible with gravity measurements of |Karatekin et aL (2006) 
(~ 6 x 10 15 kg), GRS estimation fromjKelly et al.| ( |200~6~] ) (6 x 10 15 kg) and HEND 
analysis provided by Litvak et al. (2007) ((6.1 ± 0.6) x 10 15 kg). Since our cal- 
culations are done in the domain 50°S to 90°S only and since we found the same 
numbers than other methods based on a global approach, it suggests that the CO2 
mass condensed at latitude lower than 50S is negligible. 



The GCM by Forget et al. (1999), calibrated by the pressure variation measurement 
from Viking, give also an estimate of the total SSPC of 6.9 x 10 15 kg, slightly higher 
than our estimation. The pattern of the total sublimated mass M along time (fig 
[ij is compatible with the GCM (Forget et al. (1999)) although the starting of 
sublimation seems to be slightly earlier for the D-frost model. Nevertheless, there 
is a good agreement about the end of sublimation: around Ls=270°. The shift 
in the beginning of sublimation can be explained by the integration method. In 
this study, we integrate spatially the C0 2 mass that is subliming, independently 



of the condensation domain. In the GCM integrations (Forget et al. (1999)), the 



whole SSPC spatial domain is covered: subliming and condensing CO2. At Ls=90°, 
the average SSPC balance measured by the GCM is in accumulation due to the 
geographic polar region. Our integration shows already the sublimation that occurs 
at low latitude. 



It is not possible to extend our calculation to the condensation domain because one 
energy source is missing: heat conducted in the soil (Haberle et al. (2008)), which 
is no more negligible for the accumulation period (Schmidt et al. (2009c)). Our 
estimation is limited in the domain latitude 50°S to 90°S because for lower latitudes, 
CO2 frost may be possible only with an apparent albedo below our threshold (0.3). 
Nevertheless, MOC observations shows seasonal frost at lower latitude suspected 
to be CO2 ice (Schorghofer and Edgett (2006)). These condensates could be also 
predicted with D-frost using lower albedo and the shadowing effect at lower scale 
and would be an interesting point for the future. The good agreement between our 
mass balance and direct mass measurements suggests that the low latitude frost 
mass is negligible at first order. 
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Figure 4. Total cumulative sublimated mass estimated by the D-frost model using 



OMEGA inputs compared with GCM (Forget et al.l (1999)). Please note that a neg- 



ative balance corresponds to surface mass loss (surface convention). 

If the Ar is well mixed in the south polar region, the effect on the temperature 
and thus on the total sublimated mass is below 1%. 



Figure [5] shows the integration of the SSPC sublimation for the two regions. The 
main conclusion is that despite the asymmetry of crocus lines and albedo, the 
sublimated mass seems to be symmetric around the geographic pole. The difference 
reach 20.5 % (12.5 % for t en d = U n , 31.2 % for t en d = t out ). To compile all errors, 
we estimate roughly that the mass difference is 22% ± 9. This conclusion confirms 
that the condensation phases are not equivalent during the polar night for both 



cryptic and anti-cryptic sectors in contrary of our previous work (Schmidt et al. 



(2009c)) argued only on four latitude points. There is a significant higher mass 



condensed in the anti-cryptic region, probably due to snow event, as previously 



suggested Colaprete et al. (2005); Giuranna et al. (2008). 



Assuming a density of 920 kg.m (as determined by Smith et al. (2001)), the 
estimated maximum height from the D-frost model reach 0.32 m (0.21 m for t en d = 
ti n , 0.40 m for t en d = t out ) in the cryptic region and 0.41 m (0.31 m for t en d = ti n , 



0.46 m for t en d = t out ) in the anti-cryptic region (see fig. |5J). We compile the errors 
to estimate the height 0.31 ± 0.10 m for the cryptic region and 0.38 ± 0.08 m for 
the anti-cryptic region. In comparison, the direct MOLA height measurements are 
0.37±0.08 m for the cryptic region and 0.35 ±0.08 m for the anti-cryptic region at 
Ls=150°. At first order, absolute results from both methods are compatible with a 
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difference below the error bar (see fig. [5]). The difference anti-cryptic/cryptic in the 
MOLA data is significant at Ls=210° and Ls=240° in agreement with the D-frost 
model, when the mass ratio anti-cryptic/cryptic is the highest. At this time period, 
the MOLA height measurement is higher than the D-frost estimation probably due 
to a CO2 bulk density lower than 920 kg.m -3 as suggested by Aharonson et al. 
( [2004| ; puranna et al.| ( |2007bP ; patsuo and Heki| ( |2009| >. 



During the first phase, the sublimation is nearly symmetric with a significant 
higher sublimation in the cryptic region for the time period Ls=180° to 220° (see 
fig. [6j^a) and (b )). This corresponds to appearance of the low albedo cryptic region 
( |Kieffer et aL] Q2000[ ); |Piqueux et aT| Q2003[ ); |Langevin et~aL] Q2007P ) that absorbs 
and sublimates faster. 



After Ls=220° the average anti-cryptic sector is subliming stronger with a max- 
imum difference of 1.0 x 10 8 kg.s -1 (see fig. |6^a) and (b)) despite a maximum 
sublimation rate locally in the cryptic region. This is due to the fact that the ef- 
fect of difference in surface cover is stronger than the effect of difference in albedo. 
At Ls=230, the ratio surface ratio is R surf = Surface anti- cryptic = L55 , fi u 

' "ui j bur face cryptic V L 

(a)) and the absorbance ratio is R abs = ^tlTotypT^ = °- 71 ( see fi S- @ ( b ))- 
If we simplify eq. [2] and only keep the direct sunlight term W™ n , the sublimation 
mass rates ratio can be estimated R mass = Ma M a a ^ % ~yltic %C = Rsurf -Rabs = 111, in 
agreement with our results (see fig. [6]^a) and (b)). With this simple calculation, 
we demonstrate that after Ls=220° the anti-cryptic region is sublimating faster 
due to the fact that the SSPC recession is asymmetric, i.e.: the anti-cryptic region 
reaches the outer crocus line before the cryptic region. 



Two phases are not balanced such that the total mass is different for both regions, 
with 20% more mass in the anti-cryptic region, arguing for more snow precipitation 



during winter time (Colaprete et al. (2005); Giurannaet al. (2008)). If the resulting 



material in the anti-cryptic region, more concentrated in snow, has significant lower 
emissivity, it could also condensate slower. If this process occurs during the polar 
night, 20 % is the minimum snow content. During winter time, "cold spots" has 
been observed since the 70's and been interpreted as clouds or fresh snow at the 
surface QKieffer et al] Ql977fl ; |Forget| ( |1998[ ) ; |Ivanov and Muhleman| Q2001fl ) . Recent 
studies showed that the "cold spots" dynamic in thermal infrared range is as short 



as 5 Julian days (Cornwall and Titus (2009); Cornwall and N. (2010)) suggesting 



that metamorphism is very active and very fast on Mars. Thus the CO2 emissivity 
is not likely to be lower than 0.9 during a significant period of time. 



In conclusion, our estimates indicate that the sublimation flux is modulated by 
the local properties of the cryptic/anti-cryptic regions and in particular the crocus 
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(b) Integrated mass sublimation 

Figure 5. Total sublimated CO2 mass (a) and total sublimated mass CO2 height (b) 
for both cryptic and anti-cryptic region (surface convention) assuming a density of 920 
kg.m~ 3 (as determined by Smith et al. (2001)) compared to MOLA height change. 
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(a) Absolute sublimation rates (b) Sublimation rates difference 

Figure 6. Sublimated Mass and within the cryptic region and the anti-cryptic region. Note 
that the sublimated mass is a negative quantity so that when this difference is positive, 
the cryptic region is subliming faster than the anti-cryptic region (surface convention). 



date. Schmidt et al. (2009c) conclusions are still valid: albedo is controlling the 
SSPC sublimation but once there is no CO2 ice anymore it is obvious that the 
sublimation stops. 



4 Ar dilution in the south polar vortex 



Since SSPC is sublimating, the atmosphere is pumped up and rebalanced by mass 



transfer at global scale (Hourdin et al. (1995)). This mass transfer lead also to a 
differential dilution of the minor non- condensable gas that acts like tracer of the 
phenomena. In addition to the sublimating SSPC, air mass coming from lower 



latitudes could enter the polar vortex and dilute the local air (Colaprete et al. 



(2008)). Using a box modeling, we propose here to estimate the relevant quantity 



for this system: T the advection flux and compare it to Q the sublimation flux. 



Modeling the geophysical system by boxes is a very convenient approach to under- 



stand the natural phenomena in a simple manner (Albarede (2003)). We propose 



to model the martian atmosphere using 3 boxes: the south polar deposits (SSPC), 
the polar atmosphere and the martian atmosphere as a whole (see fig [7]). 
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mars 



South polar 
atmosphere 





in 



Figure 7. Schema of the box modeling with T the advection flux, Q the sublimation flux, 
(jp i the mass of the polar atmosphere column and C the mixing ratio of non condensable 
gas, such Ar. 

4.0.0.1 C0 2 balance We could write the mass balance of the polar atmo- 
sphere of C0 2 as: 



dt 



Qin ~\~ 3~i n J~ c 



out 



(8) 



With a pol = the mass of the polar atmosphere column by surface unit (in 

kg.m -2 ), the incoming flux from the sublimation Qi n = dM estimated in the last 
section in atmosphere sign convention (in kg.m^.s^ 1 ), Ti n and Tout are respect- 
ively the incoming and outcoming advection fluxes (in kg.m _2 .s _1 ). 

We assume that the system is in steady state = such that the incoming 
fluxes exactly compensate the outgoing flux: Qi n + Ti n = J 1 out- In other terms, we 
suppose that the mass rebalance time of the atmosphere is much lower than the 
typical time of the non-condensable gas balance. 



4.0.0.2 Non-condensable species balance We could also write the mass 
balance of the minor non condensable species, assuming that they are in homo- 
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jeneous mixing within all boxes: 



With Cp i and C mars , the mixing ratio of the non condensable gas in the polar 
region respectively in the whole Martian atmosphere. We solve this equation using 
the following numerical solution: 

C pol (t + At) = C pol (t). (l - —.At - ^.At) + C mars .—.At (10) 



For the convenience, we use a time step At = Is. 

Measurements of the Ar in the south polar atmosphere have been done by the GRS 



instrument (Sprague et al. (2004); Sprague et al. (2007)). During the time period of 



sublimation in the south polar region (Ls= 11 0-350), we estimate the decrease of Ar 
using eq. 10 We set the initial value of Ar abundance C po i(t = 110°Ls) = 3.89% as 



measured by GRS. We use the commonly used value of the Martian mixing ratio 



Cmars = 1-6% in agreement with Viking Lander 2 measurements (Owen et al. 



(1977)). 



Minor species abundances evolution with time are plotted in fig. |SJ for different 
values of Ti n : 0, 4.10 -6 , 8.10 -6 and 20. 10~ 6 kg.m _2 .s _1 . The mixing ratio evolution 
decreases with time to reach a minimum around Ls=250 and then increases again. 
The first decrease is due to the sublimation of the seasonal frost, free of non- 
condensable gas. In addition, the advection flux dilutes Ar with a higher Ti n leading 
to a faster dilution. The increase in second step is due to the injection of the 
Martian air, slightly enriched in Ar compared to the local polar air. The rate of 
increase is higher when Ti n is faster. 

By matching GRS data with our box model, we estimate that the flux J 7 ^ ~ 
4.10 -6 kg.m _2 .s _1 . This value has the same order of magnitude as the mean value 
of the sublimation flux (Qi n ) = 6.7 x 10~ 6 kg.m~ 2 .s _1 . This result suggests that 
the air mass entering the south polar vortex during the recession phase is of the 
same order than the sublimated mass. Thus the local difference in Qi n due to 
asymmetry should lead to non-negligible effect on the vortex dynamic, especially 
after Ls=220° (see fig.g). 

Future studies using GCM or local climate model have to be conducted to study the 
actual dynamical effect of the asymmetry in sublimation. One particular attention 
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Figure 8. Estimation of C po i as a function of time in comparison with actual Ar measure- 
ments in the South Polar atmosphere. Lines represent the solution of our box modeling 
(eq. 10) for four cases: Ti n = 0, 1CP 5 , 5.10 -5 and 10 -4 kg.m~ 2 .s . Triangles are Ar 
observations of the GRS instruments (Sprague et al. ( 2007| )) for the whole south polar 
region. 

should be drawn on the dilution factor regarding the dynamics in the cryptic vs 
anti-cryptic region. Some regional differences should be pointed out, as suggested 
by GRS (|Sprague et aT] (|2007j)). 



5 Conclusions 



We estimate the CO2 mass balance of the SSPC using the D-frost model and 
estimate the total sublimated mass to be around (5.6 ±1,3) x 10 15 kg which is 



compatible with gravity measurements (Karatekin et al. (2006)), with measure 



ments by gamma ray and neutron spectroscopy (Litvak et al. (2006); Kelly et al. 



(2006)) and with previous GCM studies (Forget et al. (1999); Kelly et al. (2006)). 



This agreement validates our current approach. 

The effect of the Ar content, in the case of a well mixed atmosphere is below 1% 
on the total sublimated mass. A more realistic case of denser Ar enriched air in 
the bottom part of the atmosphere should be studied in the future. 



We show that the SSPC is not symmetric in mass around the geographic pole, i.e.: 
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the cryptic and anti-cryptic regions have an accumulated mass with a difference 
in order of 22 ± 9 %. in favour of the anti-cryptic region. Previous studies sugges- 
ted that in addition to direct condensation, during the southern fall and winter 
seasons, snow precipitation is more intense in the anti-cryptic region than in the 



cryptic one (Colaprete and Toon (2002); Giuranna et al. (2008)). Our results show 
that, the actual mass deposited (possibly as snow) should be in order of 20 %. 
This conclusion is still valid if the CO2 emissivity has regional and/or temporal 
variations within the range 0.9 to 1 during the spring and summer, as observed by 



TES ( Eluszkiewicz and Moncet (2003)) ans supported by fast CO2 metamorphism 



( Cornwall and N. ( 2010)). T he previous estimation on small regions of interests 
(ROI) by|Kieffer et al.| 02000} was 852 kg.m" 1 in the "cryptic" ROI and 841kg.ni" 1 



in the "bright cap" ROI. The difference is probably due to the crocus date meas- 
urement that is delayed for OMEGA in comparison to TES. In case of subpixel 
mixing, due to the thermal emission dependence on the temperature at power four, 
bare soil is dominating in the thermal infrared. This effect is not present in near 
infrared domain and the mixing is really linear. Further studies should be done to 
precisely compare visible, near infrared and thermal infrared datasets. 

From Ls=180° to 220°, we point out that the cryptic region is sublimating slightly 
stronger than the anti-cryptic region. In the second phase, anti-cryptic region 
sublimation is dominating with a maximum difference in sublimation of 1.0 x 10 8 
kg.s" 1 . This results is compatible with the MOLA height measurements difference 
that reach a maximum for Ls=210°-240°. 

We use a simple box model to fit the GRS measurements of Ar and estimate that 
the atmosphere advection flux entering the south polar vortex is similar to the flux 
from the sublimating SSPC. 



More precise studies on atmosphere dynamics and minor species concentration 
should be carried out with a GCM or a local scale climate model. At present time, 
SSPC albedo is a constant parameter in any GCM. Future studies should be done 
to include the SSPC asymmetry of albedo leading to an asymmetry of sublimation 
in order to simulate realistic dynamics and dilution factors. Those results will have 
to be compared to the atmospheric observations of pressure (for instance using the 



(using gamma rays, i.e.: 



method proposed by Forget et al. (2007); Spiga et al. (2007)) and minor species 
concentration 



on the ground i.e: Encrenaz et al. (2005), or orbiters i.e.: Melchiorri et al. (2007)). 



Sprague et al. (2004), optical measurements 



In particular, future local climate studies should explain the actually asymmetric 



Ar enrichment (Sprague et al. (2004)) that could play an indirect role in the 



stability of the SSPC itself. If the atmosphere is locally enriched in Ar, the partial 
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pressure of CO2 is lower and then CO2 ice stability is lower (Forget et al. 



(2008)). 
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